(C--,  ?■' 


vi 


r  j,i  U  , 

t  i 


p'-ll 

IMI 


1  e  ^  ^ 


SLL  80  619 

copy  1 

,  ,X^  ^7(i,  f 


Lr^e.!^ 


MODELING  C.W.  DF  AND  HF  LASER  PERFORMANCE 
AT  LOW  CAVITY  PRESSURES  AND  DILUENT  CONDITIONS 

W.L.  Rushmore  and  S.W.  Zelazny 


NIAGARA  FRONTIER  OPERATIONS 


Bell  Aerospace 


TEXTRON 


Division  of  Textron  Inc. 


Buffalo,  New  York  14240 


AIAA  Paper  No.  IV-5 

AIAA  Conference  On 
Fluid  Dynamics  of  High  Power  Lasers 
Cambridge,  Massachusetts 

OCTOBER  31  -  NOVEMBER  2,  1978 


PLEASE  RETURN  TO: 

BMD  TECHNICAL  INFORMATION  CENTER 
BALLISTIC  MISSILE  DEFENSE  ORGANIZATION 
7100  DEFENSE  PENTAGON 
WASHINGTON  D.C.  20301-7100 

-  '--s  ^ 


^  Approved  for  public  release; 
,  Distribution  Unlimited 


■O.HoT-'I' 


Accession  Number:  4022 


Publication  Date:  Nov  02, 1978 

Title:  Modeling  CW  DF  and  HF  Laser  Performance  at  Low  Cavity  Pressures  and  Diluent  Conditions 
Personal  Author:  Rushmore,  W.L.;  Zelazny,  S.W. 

Corporate  Author  Or  Publisher:  Bell  Aerospace  TEXTRON,  Buffalo,  NY  14240  Report  Number:  AlAA 
Paper  No.  IV-5  Report  Number  Assigned  by  Contract  Monitor:  SLL  80  619 

Comments  on  Document:  Archive,  RRl,  DEW.  Presented  at  the  AIAA  Conference  on  Fluid  Dynamics  of 
High  Power  Lasers,  Cambridge,  MA,  October  31 -November  2, 1978 

Descriptors,  Keywords:  Model  Continuous  Wave  Hydrogen  Deuterium  Huoride  Laser  Performance  Low 
Cavity  Pressure  Diluent  Condition  Geometry  Flow  Area  Expansion  Chemiluminescence  Vacuum  Duct 
Pressure  Chemical 

Pages:  00039 

Cataloged  Date:  Dec  07,  1992 
Document  Type:  HC 
Number  of  Copies  In  Library:  000001 
Record  ID:  25501 


Source  of  Document:  DEW 


MODELING  C.W.  DF  AND  HF  LASER  PERFORMANCE 
AT  LOW  CAVITY  PRESSURES  AND  DILUENT  CONDITIONS* 


W.L.  Rushmore^  and  S.W.  Zelazny 

BalABraspanDlilllSI!]  Buffalo,  New  York 

Division  of  Textron  Inc. 

Abstract 

A  laser  performance  model  originally  developed  for  high  cavity  pressure  and  diluent  c.w. 
lasers  is  extended  to  apply  to  a  significantly  lower  pressure  and  diluent  operating  regime.  The  dif¬ 
ferent  class  of  laser  nozzle  geometries  used  in  low  diluent  applications  are  described  and  a  method 
for  theoretically  modeling  the  internal  flow  area  expansion  which  occurs  in  these  designs  is  presented. 
Experimental  data  (closed  cavity  power,  chemiluminescence,  zero  power  gain)  for  helium  diluent 
hydrogen  fluoride  and  nitrogen  diluent  deuterium  fluoride  chemical  lasers  is  used  to  assess  the 
accuracy  of  the  model.  Agreement  between  theory  and  data  is  good  for  both  the  axisymmetric  and 
two  dimensional  nozzle  geometries  examined.  Test  conditions  where  theory  and  data  do  not  agree 
are  shown  to  be  in  part  related  to  the  influence  of  shroud  geometry  and  vacuum  duct  pressure  on 
the  flowfield. 

Nomenclature 

A  =  Laser  nozzle  face  area 

m  =  Mass  flow 

,  *  • 

^^F2^A  “  Total  molar  flowrate  of  available  fluorine,  Np/2  +  Np^ 

Nj  =  Molar  flowrate  of  specie  i 

p7  =  Recoverable  pressure  at  diffuser  exit 

Presented  at  the  AIAA  Conference  on  Fluid  Dynamics  of  High  Power  Lasers,  Cambridge,  Mass., 
October  31  ,-November  2,  1978. 

Index  Categories;  Lasers:  Reactive  Flows 
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=  Combustor  pressure 

=  Nozzle  exit  pressure 

=  Total  outcoupled  power  for  DF*  lasing 

=  Outcoupled  power  for  lasing  from  v  to  v-1  vibrational  level 

=  Small  signal  gain  for  lasing  from  v+1  to  v  vibrational  level,  on  rotational  line  J 

=  Reynolds  number  based  on  stagnation  conditions  and  nozzle  radius  or  half  height 

=  Cavity  molar  mixture  ratio 

=  Combustor  total  temperature 

=  Total  temperature  at  injector  exit 

=  Axial  velocity 

=  Axial  coordinate 

=  Axial  location  of  optimum  stable  resonator  axis 
=  Combustor  molar  diluent  ratio 
=  Combustor  molar  diluent  ratio  of  N2 
=  Cavity  molar  diluent  ratio 
=  Power  flux  per  unit  area 
=  Total  (combustor  plus  cavity)  diluent  ratio 

=  Specific  power  based  on  combustor  and  lasant  nozzle  mass  flow 
=  Molar  specific  power/reference  molar  specific  power,  see  ref.  2  for  reference  value. 


I.  Introduction 

The  Problem 

The  basic  components  of  a  chemical  laser  are  shown  schematically  in  Figure  1 .  In  an  earlier 
study  a  theoretical  model  was  described  which  predicted  the  closed  cavity  power  in  continuous 
wave  (c.w.)  HF  and  DF  lasers.  The  operating  regime  examined  in  that  study^  was  generally  charac¬ 
terized  by  mass  throughputs,  in/A,  greater  than  1 .0  g/s-cm^  and  high  diluent  levels  >  30).  Also, 
the  laser  nozzle  designs  considered  in  this  earlier  study  ^  contained  very  little  internal  area  relief  for 
flow  expansion  induced  by  heat  release  in  the  laser  cavity.  Comparison  of  the  schematics  of  the  low 
base  relief  BCL-7  and  high  base  relief  BCL-1 3  nozzles  shown  in  Figure  2  illustrates  the  difference 
between  the  two  types  of  nozzle  concepts. 
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The  high  base  relief  nozzle  has  practical  application  in  laser  systems  which  must  achieve 
high  power  extraction  efficiency  at  low  diluent  levels.  In  addition,  these  laser  systems  only  require 
relatively  low  pressure  recovery  capability  (30  torr  or  less).  The  specific  problem  addressed  in  this 
study  was  to  evaluate  the  capability  of  the  laser  performance  analysis  technique*  developed  for 
high  pressure,  high  diluent  lasers  to  characterize  this  significantly  different  operating  regime  and  to 
improve  the  model  where  necessary. 

Approach 

A  number  of  investigations'^'"  have  been  made  which  provide  closed  cavity  power,  chemilum¬ 
inescence  measurements,  and  zero  power  gain  data  for  the  operating  regime  of  interest,  i.e.,  low 
diluent  and  cavity  pressures.  Two  approaches  have  been  used  to  evaluate  the  ability  of  a  laser  model 
to  characterize  details  of  the  laser  flowfield.  The  first  approach  uses  one  particular  nozzle  geometry 
and  operating  condition ',  whereas  the  second  approach  uses  a  number  of  nozzle  geometries  and 
operating  conditions*.  We  have  selected  the  second  approach  in  this  study.  The  governing  equations, 
method  of  solution,  and  technique  used  to  establish  the  conditions  at  the  cavity  entrance  plane  are 
identical  to  those  described  by  Zelazny  et  al  * . 

The  data  base'^'"  used  in  this  study  is  described  in  Section  11.  Details  of  the  analyses  used 
to  consider  influences  of  the  base  area,  stable  resonator  geometry  and  cavity  pressure  distributions 
are  given  in  Section  III.  Results  are  given  in  Section  IV  followed  by  the  main  conclusions  of  this 
investigation.  Section  V. 

II.  Data  Base 

A  summary  of  the  data  reviewed  in  this  study  is  given  in  Table  1 .  The  axisymmetric  nozzles 
BCL-13  and  BCL-14  are  shown  schematically  in  Figure  3,  whereas  the  two-dimensional  nozzle  LRC-1 
and  quasi-axisymmetric  nozzle  LRCL-2  are  shown  in  Figure  4.  The  LRCL-2B  nozzle  differs  from 
the  LRCL-2  nozzle  in  that  the  sonic  lasant  injection  orifices  of  LRCL-2  were  modified  to  supersonic 
conical  nozzles  in  the  LRCL-2B  configuration^.  This  design  (LRCL-2  and  2B)  is  herein  referred  to 
as  a  quasi-axisymmetric  nozzle  since  unlike  the  BCL-13  and  14  designs  it  does  have  azimuthal  vari¬ 
ations  in  the  lasant  flowfield  around  the  circumference  of  the  circular  primary  nozzle  due  to  the 
discrete  orifice  injection  configuration. 
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The  data  to  be  examined  can  be  conventiently  subdivided  into  helium  diluent  hydrogen 
fluoride  data  and  nitrogen  diluent  deuterium  fluoride  data  designated  as  He-HF*  and  N2-DF*. 

Three  types  of  data  are  available  (1)  closed  cavity  power  (CCP)  obtained  with  a  stable  resonator 
which  is  used  to  determine  the  optimum  optical  axis  (x^)  and  provide  an  approximate  measure  of 
the  lasing  zone  length  (2Xp);  (2)  chemiluminescence  data  giving  the  number  density  and  rotational 
temperature  distribution  in  the  lasing  zone  and  (3)  zero  power  gain  distributions  in  the  flow  direction 
on  the  nozzle  centerline. 

The  laser  performance  may  be  expressed  in  terms  of  the  following  flow  parameters;  (1)  mass 
throughput,  m/A,  (2)  combustor  diluent  level,  (3)  adiabatic  flame  temperature  in  the  combustor, 
TqJ,  (4)  cavity  mixture  ratio,  R^,  and  (5)  diluent  level  in  the  lasant  stream,  Examination  of  the 
experimental  data  base  shows  how  performance  varies  with  these  parameters  for  a  given  nozzle  geo¬ 
metry.  Results  of  the  theoretical  model  predictions  are  compared  with  the  experimentally  observed 
trends  in  Section  IV. 

III.  Analysis 

Either  two-dimensional  or  axisymmetric  flowfields  are  considered  and  the  boundary  layer 
approximations  are  assumed  to  apply.  A  quasi-two-dimensional  form  of  these  equations  were 
obtained^  and  solved  using  implicit  integration  techniques  by  the  BLAZE-II  Computer  Code. 

Starting  conditions  at  the  nozzle  exit  plane  were  computed  using  the  CNCDE  computer  code 

O 

developed  by  Driscoll  where  viscous  losses  induced  by  the  nozzle  boundary  layers  are  considered 
for  both  the  primary  and  secondary  nozzles.  Combustor  heat  loss  and  totally  regenera tively  cooled 
nozzle  designs,  e.g.,  BCL-13,  are  considered  by  transferring  the  primary  nozzle  heat  loss  to  the 
secondary  nozzle  to  raise  the  lasant  stagnation  temperature.  The  following  subsections  describe 
features  of  the  analysis  which  have  been  included  to  account  for  nozzle  designs  used  in  the  low 
pressure  and  diluent  operating  range  of  interest  in  this  investigation. 

Considering  Base  Area  Effects 

Figure  5  shows  a  schematic  of  the  side  view  of  an  axisymmetric  or  two-dimensional  (slit) 
nozzle  which  contains  a  significant  base  region.  A  number  of  studies,  e.g.  ref  9,  have  examined  the 
details  of  the  flowfield  which  include  the  influence  of  transverse  pressure  gradients,  flow  separation 
of  the  primary  and/or  secondary  nozzle  and  the  recirculation  zone  in  the  base  region.  Unfortunately, 
the  elliptic  equation  system  required  to  include  these  influences  requires  considerable  computer 
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solution  time  and  the  resulting  models  cannot  be  readily  used  to  conduct  extensive  parametric 
studies  to  establish  the  trends  needed  to  guide  experimental  test  plan  definition.  The  model  used 
herein  uses  a  parabolic  equation  system  and  assumes  the  pressure  across  the  mixing,  reacting,  and 
lasing  region  is  constant.  Hence,  transverse  pressure  gradients  and  recirculation  effects  have  been 
neglected.  As  a  consequence  of  neglecting  these  effects,  it  becomes  necessary  to  develop  a  technique 
which  allows  a  reasonable  definition  of  effective  (or  idealized)  nozzle  exit  plane  conditions.  Speci¬ 
fically,  the  flow  and  geometric  parameters  which  must  be  defined  are  an  average  static  pressure, 
effective  primary  and  secondary  nozzle  exit  dimensions  (the  mixing  scales)  and  temperature  and 
velocity  boundary  layer  profiles  which  reflect  the  flow  expansions  and/or  compressions  consistent 
with  the  average  pressure.  These  parameters  are  computed  as  follows.  For  a  sonically  injected 
lasant  stream,  the  flow  is  assumed  to  expand  outwards  into  the  base  region  until  its  pressure  is 
matched  with  the  primary  nozzle  exit  pressure.  This  implies  that  the  primary  nozzle  exit  pressure 
is  the  critical  pressure  for  the  cavity  mixing  calculation.  The  initial  laminar  mixing  length  scale  for 
the  secondary  nozzle  is  then  taken  to  be  the  sum  of  the  expanded  secondary  nozzle  dimension  plus 
the  remaining  base  dimension.  As  will  be  shown  in  Section  IV,  these  assumptions  allow  the  theo¬ 
retical  predictions  for  the  main  flow  features  and  gain  distributions  to  be  accurately  correlated  with 
the  experimental  observations.  It  is  the  treatment  of  the  internal  base  region  which  accounts  for 
the  principal  difference  in  modeling  high  pressure  (as  in  ref.  1 )  and  low  pressure  flows.  As  the  inter¬ 
nal  base  approaches  zero,  the  two  models  become  equivalent. 

The  pressure  distribution  in  the  laser  cavity  is  a  function  of  the  imposed  boundary  condition. 
Unless  specifically  identified,  the  results  presented  herein  assume  a  constant  pressure  cavity  until 
the  internal  base  is  filled.  From  this  point  on,  a  constant  area  flow  (equal  to  the  laser  nozzle  face 
area  A)  is  assum^  and  hence  a  subsequent  pressure  rise  is  computed.  It  should  be  pointed  out  that 
a  freejet  flow,  i.e.,  an  unshrouded  laser  flow  which  surrounds  the  laser  pressure  field  with  the  vacuum 
duct  pressure,  is  not  equivalent  to  a  constant  pressure  flow  in  that  the  bulk  average  pressure  in  a 
freejet  flow  can  decrease  if  the  exit  flow  is  underexpanded. 

Closed  Cavity  Power 

The  closed  cavity  power  is  computed  using  a  Fabry-Perot  optical  cavity  and  the  gain  equals 
loss  assumption  is  made.  The  mirror  reflectivities  are  specified  as  input  parameters.  A  laminar 
mixing  model  is  used  which  takes  into  account  the  local  pressure  and  temperature  to  compute  a 
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laminar  diffusion  coefficient.  Rotational  equilibrium  is  assumed  for  all  vibrational  states  (see 
Sentman^®  for  a  discussion  on  the  consequences  of  this  assumption).  The  DF  kinetics  are  the 
same  as  those  used  in  ref.  1  and  are  given  in  Table  2.  The  HF  kinetics  are  those  given  in  ref.  11, 
and  are  listed  in  Table  3.  Calculations  were  made  assuming  either  full  F-atom  recombination  or 
negUgible  F-atom  recombination  on  the  nozzle  wall  at  the  cavity  entrance  plane.  The  effect  of  the 
F-atom  recombination  assumption  will  be  shown  in  the  results  section.  Upstream/downstream 
coupling  effects  induced  by  the  stable  resonator  configuration  which  reduce  laser  closed  cavity  power 
due  to  the  propagation  of  radiation  through  regions  of  negative  gain  (absorption  zones)  are  obviously 
not  considered  in  the  Fabry-Perot  model.  An  approximate  method  of  accounting  for  this  influence 
in  this  model  will  be  described  in  the  following  section. 

rV.  Results 

N2-DF*  Data  and  Theory  Comparisons 

Closed  cavity  power  and  chemiluminescence  data  was  described  by  Zelazny  et  al  for  the 
two  axisymmetric  nozzles  BCL-13  and  14  shown  in  Figure  3.  Both  nozzles  have  identical  primary 
and  secondary  dimensions  and  differ  only  in  the  nozzle  packing  density.  The  BCL-13  nozzle  has 
45  nozzles/in^,  whereas  BCL-14  has  33  nozzles/in^.  Consequently,  the  BCL-14  nozzle  provides 
12%  more  base  region  for  flow  expansion  than  the  BCL-13  array  but  requires  combustor  pressures 
36%  greater  than  BCL-14  to  achieve  the  same  mass  flowrate. 

Figure  6  shows  the  rotational  temperature  distribution  with  distance  downstream  obtained 
from  the  chemiluminescence  data  for  BCL-13.  Also  shown  are  the  predictions  using  the  C/B  (i.e., 
CNCDE-BLAZE-Il)  computer  code  with  various  methods  considered  for  treating  the  internal  base. 
The  baseline  curve  is  obtained  using  the  assumptions  outlined  in  Section  III,  i.e.,  the  secondary 
flow  expands  to  a  matched  pressure  and  it’s  laminar  mixing  scale  length  is  the  entire  base  dimension. 
Two  other  methods  for  modeling  the  base  region  were  examined,  and  results  are  also  shown  in 
Figure  6.  In  one  method  the  laminar  mixing  length  scale  of  the  secondary  nozzle  was  taken  as  that 
dimension  required  to  obtain  a  matched  pressure  with  the  primary  nozzle  flow.  The  second  method 
used  the  entire  base  region  for  complete  expansion  of  the  nozzle  flow  at  the  nozzle  exit  plane  before 
the  reaction  commences.  A  matched  pressure  plus  primary  and  secondary  length  scales  are  then 
determined  based  on  the  total  nozzle  bank  area  available  to  the  flow.  As  can  be  seen,  the  baseline 
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curve  follows  the  data  much  more  closely,  particularly  in  the  first  centimeter  where  the  majority 
(>  75%)  of  the  power  is  extracted. 

Detailed  closed  cavity  power  data  for  BCL-14  was  used  to  compare  theory  with  data.  Para¬ 
metric  studies  were  conducted  using  the  BLC-14  nozzle  varying  p7  (i.e.,  the  recoverable  pressure 
assuming  80%  normal  shock  recovery  in  a  constant  area  supersonic  diffuser),  combustor  inlet 
temperature  Tqj,  combustor  molar  diluent  level  ,  cavity  molar  mixture  ratio  Rl,  and  cavity  molar 
diluent  ratio  |3g.  The  recoverable  pressure  is  directly  correlatable  with  the  mass  throughput  at  a 
fixed  TqJ,  Rj^,  and  as  shown  in  Figure  7  for  the  BCL-13  and  14  nozzle  designs. 

Figure  8  shows  the  predicted  power-on  gain  profiles  for  BCL-14.  The  gain  curves  are  for 
transitions  from  the  v=l  to  v=0  lines  for  various  J  values.  As  a  consequence  of  the  rotational 
equilibrium  assumption,  the  power  is  computed  such  that  extraction  can  only  occur  from  one  J 
value  at  a  time  for  each  vibration  line^^.  As  a  result,  the  lower  J  values  enter  an  absorption  condi¬ 
tion  while  power  is  being  extracted  from  the  higher  J  levels. 

This  neglecting  of  multiline  lasing  and  the  potential  for  the  radiation  intensity  propogated 
from  a  positive  gain  region  through  an  absorption  zone  downstream  represents  a  significant  short¬ 
coming  of  the  current  model.  The  result  of  this  limitation  is  that  laser  power  predictions  made 
assuming  an  optical  cavity  composed  of  plane  parallel  mirrors  will  allow  power  extraction  on  lines 
which  will  not  achieve  saturation  in  stable  resonators  containing  concave  mirror  geometries.  A 
very  simple  yet  effective  method  for  approximating  this  effect  without  resorting  to  costly  (increased 
computer  time)  modifications  to  the  resonator  model  is  to  reduce  the  power  extraction  efficiency 
of  the  Fabry-Perot  cavity  prediction.  Here,  we  uniformly  downgrade  the  power  extraction  obtained 
on  the  v=3-2,  2-1,  and  1-0  transitions  by  15%  based  on  comparisons  between  theory  and  data. 

Figure  9  shows  the  C/B  comparisons  with  BCL-14  data  for  a  P7  (i.e.,  m/A)  scan.  Two  pre¬ 
diction  curves  are  presented.  One  was  made  assuming  complete  F-atom  recombination  on  the 
nozzle  wall,  whereas  the  second  was  made  assuming  no  F-atom  recombination.  As  can  be  seen. 


the  complete  recombination  results  agree  much  better  with  the  F*  and  x^  data.  For  this  reason, 
the  remainder  of  the  C/B  N2-DF*  predictions  are  made  with  complete  recombination.  Figure  10 
shows  the  results  for  a  T^j  scan.  The  predictions  correctly  show  the  increase  of  F*  with  T^j  to 


within  10%  of  the  data.  Figure  1 1  shows  the 


scan  results.  The  data  shows  that  a*  remains 
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essentially  constant  as  is  raised  from  1.0  to  8.0,  while  the  C/B  predictions  show  a  20%  drop  in 
'a*.  However,  the  predictions  do  correctly  show  the  increase  in  with  Figure  1 2  shows  the 
results  of  the  scan.  BCL-14  was  designed  to  operate  in  either  the  regeneratively  cooled  mode 
(i.e.,  the  secondary  nozzle  flow  is  used  to  cool  the  complete  nozzle  bank)  or  the  dump  cooled  mode 
(where  some  of  the  secondary  flow  is  dumped  externally  from  the  cavity).  In  this  way  low  Rl 
values  could  be  run  in  the  dump  cooled  mode.  For  the  data  of  Figure  12,  the  Rj^  =  5.0  and  10.0 
conditions  were  run  in  the  dump  cooled  mode.  The  CNCDE  code  may  also  be  run  in  the  dump 
cooled  mode  (i.e.,  the  heat  transfer  from  the  primary  flow  to  the  secondary  flow  is  controlled),  so 
the  C/B  predictions  for  Rj^  =5.0  and  10.0  were  also  made  in  the  dump  cooled  mode.  As  can  be 
seen,  the  increase  in  a*  as  Rj^  decreases  is  correctly  predicted  to  within  18%.  Finally,  Figure  13 
shows  the  results  for  the  scan.  The  data  and  theory  both  show  no  benefit  in  increasing  and  an 
insensitivity  of  lasing  zone  length  to 

Comparisons  between  theory  and  data  for  the  BCL-13  nozzle  showed  results  similar  to  those 
obtained  with  the  BCL-14  nozzle  comparisons.  The  theory  showed  that  the  BCL-13  design  would 
yield  approximately  15%  higher  efficiency  than  the  BCL-14  nozzle  at  comparable  throughput  con¬ 
ditions.  Table  4  shows  a  comparison  between  theory  and  data  for  the  BCL-13  nozzle  at  a  high 
performance  condition. 

He-HF*  Data  and  Theory  Comparisons 

Closed  cavity  power  data  for  the  LRCL-1 ,  2,  and  2B  nozzles  in  a  freejet  configuration  were 
examined.  The  LRCL-1  two-dimensional  and  LRCL-2  axisymmetric  nozzle  geometries  are  shown 
in  Figure  4.  The  LRCL-2  nozzle  is  different  from  the  BCL-13  and  BCL-14  axisymmetric  nozzles  in 
that  the  LRCL-2  secondary  nozzles  consist  of  six  sonic  orifices  equally  spaced  around  the  primary 
nozzle,  whereas  the  BCL-13  and  1 4  axisymmetric  secondary  nozzles  consist  of  a  concentric  circular 
orifice  (see  Figure  3).  The  LRCL-2  nozzle  was  modeled  as  a  “true”  axisymmetric  in  that  an  equiva¬ 
lent  area  concentric  secondary  nozzle  was  assumed.  The  LRCL-2B  nozzle  has  an  identical  primary 
nozzle  to  LRCL-2  with  six  supersonic  secondary  nozzles  with  an  area  ratio  of  16.0.  Again,  this  was 
modeled  as  an  equivalent  throat  and  exit  area  concentric  ring. 

Comparisons  of  C/B  predictions  with  the  zero  power  gain  data  for  LRCL-2  are  shown  in 
Figures  14  and  15  for  the  P^fb)  and  P2(5)  lines.  The  effect  of  F-atom  wall  recombination  on  the 
nozzle  wall  at  the  exit  plane  is  also  examined.  As  was  the  case  for  the  BCL-14  power  data,  the 
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complete  F-atom  recombination  predictions  follow  the  data  much  more  closely.  The  C/B  runs 
were  made  assuming  a  constant  area  cavity  boundary  condition.  The  internal  base  was  sufficient 
so  that  all  runs  predicted  a  constant  pressure  flow  to  5  cm.  The  data  was  taken  in  a  freejet  con¬ 
figuration,  and  as  stated  previously,  this  may  not  be  equivalent  to  a  constant  pressure  cavity  flow 
if  the  vacuum  duct  pressure  is  lower  than  the  nozzle  exit  pressure.  This  may  explain  why  the  Pjfb) 
data  has  a  longer  positive  gain  region  than  the  predictions.  The  P2(5)  data,  which  goes  into  absorp¬ 
tion  sooner,  is  closely  tracked  by  the  C/B  predictions. 

Figure  16  shows  the  number  density  data  for  HF(1)  and  FIF(2)  compared  with  the  C/B 
predictions.  Again,  the  complete  recombination  results  track  the  data  much  more  closely  than  the 
no  recombination  results,  particularly  for  HF(2).  Therefore,  the  remainder  of  the  C/B  runs  for 
He-HF*  are  also  made  with  complete  F-atom  wall  recombination.  Figure  1 7  shows  the  zero  power 
gain  results  for  LRCL-1  for  the  P|  (7)  and  P2(”7)  lines.  While  the  C/B  predictions  are  not  in  as  good 
agreement  as  for  LRCL-2,  the  peak  gains  are  still  predicted  to  within  1 .5%/cm  and  the  location  of  the 
peak  gain  to  within  0.5  cm. 

The  arguinent  made  for  N2-DF  lasers  that  upstream/downstream  coupling  effects  due  to 
absorption  will  degrade  the  actual  power  is  also  relevant  to  the  HF  data  examined.  Hence,  the  C/B 
power  predictions  were  degraded  by  25%  to  give  good  agreement  between  theory  and  data.  Specifi¬ 
cally,  the  following  ralation  was  used. 

pHF*=  0.75  (PHF%  pHF*) 

Tot  1.0  2-1 

It  is  seen  from  the  above  relationship  that  for  DF  lasing  the  power  available  is  obtained  from 
three  levels,  whereas  for  HF  lasing  it  is  available  from  only  two  levels. 

The  freejet  data  for  LRCL-1 ,  2  and  2B  was  examined  and  the  results  are  presented  in 
Figure  1 8.  The  LRCL-2B  model  requires  some  further  discussion.  As  stated  in  Section  HI,  the 
sonic  secondary  flow  was  allowed  to  expand  into  the  base  region  until  a  matched  pressure  with  the 
primary  nozzle  flow  is  obtained.  For  the  supersonic  secondary  LRCL-2B  nozzle,  the  secondary 
nozzle  exit  pressure  was  typically  only  10  to  20%  that  of  the  primary  nozzle.  Since  the  BLAZE-II 
code  cannot  consider  transverse  pressure  gradients  explicitly,  the  initial  pressure  for  the  laminar 
mixing  calculations  was  taken  as  the  primary  nozzle  pressure  since  the  dominant  mixing  mechanism 
is  the  diffusion  of  H2  into  the  primary  flow. 
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At  the  lower  flowrates  (m/A  =  0.05  and  0.10),  the  C/B  predictions  are  within  10%  of  the 
data  and  correctly  predict  that  LRCL-2  gives  the  best  performance  (i.e.,  a),  with  LRCL-1  and 
LRCL-2B  giving  almost  equivalent  specific  power  a.  At  the  higher  flowrates,  the  C/B  predictions 
start  falling  off  and  give  significantly  lower  a’s  than  measured.  A  possible  explanation  for  this  result 
is  that  in  a  freejet  configuration  the  influence  of  the  lower  vacuum  duct  pressure  will  have  more  of 
an  influence  at  higher  m/A  (i.e.,  higher  exit  pressures).  The  initial  C/B  runs  were  made  with  a  con¬ 
stant  area  cavity  boundary  condition.  The  LRCL-2  and  LRCL-2B  predictions  show  adequate  inter¬ 
nal  base  to  allow  constant  pressure  cavity  flow  at  all  flowrates.  However,  the  LRCL-1  predictions 
show  that  the  flow  will  choke  at  the  higher  m/A’s  of  0.29  and  0.35  g/cm^.  Therefore,  these  con¬ 
ditions  were  rerun  with  a  linear  area  expansion  of  15%  over  10  cm,  which  was  adequate  to  predict 
a  constant  pressure  cavity  flow,  see  Figure  1 8. 

The  effect  of  the  cavity  boundary  conditions  are  further  examined  in  Figure  19.  Two 
additional  sets  of  calculations  for  LRCL-2  were  made.  Both  assumed  that  the  cavity  pressure 
decreases  from  the  nozzle  exit  pressure  to  a  value  of  1  torr  at  5  cm  downstream.  One  set  of  runs 
assumed  a  linear  decrease,  while  the  second  set  assumed  a  parabolic  decrease  (i.e.,  the  pressure  drops 
off  faster  in  the  initial  region  than  the  case  using  a  linear  profile.)  As  can  be  seen,  the  predictions 
with  the  faster  (parabolic)  decrease  in  pressure  show  the  same  trend  as  the  data  in  that  the  specific 
power  no  longer  drops  off  sharply  with  increasing  m/A.  Thus  the  cavity  boundary  condition  has  a 
marked  effect  when  considering  the  scalability  of  the  data  to  larger  nozzle  dimensions.  Table  5 
shows  a  comparison  of  the  measured  and  predicted  for  the  data  of  Figures  18  and  19.  As  can 
be  seen,  the  predicted  Xj,’s  are  typically  40-70%  shorter  than  the  measured  values  when  a  constant 
area  boundary  condition  is  employed;  however,  when  the  parabolic  decrease  in  pressure  boundary 
condition  was  used,  the  predicted  x^’s  for  LRCL-2  increase  and  are  much  closer  to  the  measured 
values. 

A  number  of  uncertainties  in  the  accuracy  of  the  reaction  rate  models  exist  and  hence  it  was 
of  interest  to  determine  the  sensitivity  of  computed  laser  performance  to  changes  in  various  rate 
constants  for  a  low  pressure  He-HF  laser.  The  test  conditions  chosen  were  for  LRCL-2  at  a  low  m/A, 
since  both  the  predicted  gain  and  performance  were  both  in  good  agreement  with  experimental  data. 
The  results  are  given  in  Table  6.  Different  sets  of  rates  were  examined  individually  with  respect  to 
the  baseline  rates.  Multiplying  or  dividing  the  pumping  rates  by  a  factor  of  two  had  the  expected 
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result  of  increasing  and  decreasing  the  predicted  specific  power.  In  all  cases  doubling  the  deacti¬ 
vation  rates  also  decreased  specific  power.  However,  the  net  effect  was  always  less  than  10%. 

V.  Conclusions 

(1)  A  direct  extension  of  the  model  developed  in  low  base  nozzles  designed  for  hi^  cavity 
pressure  operation^  was  found  to  inadequately  characterize  the  rate  of  mixing  between  the  lasant 
and  primary  streams  of  high  base  nozzles.  Numerical  experimentation  showed  that  the  mixing  scale 
required  to  accurately  describe  the  entrainment  of  the  lasant  into  the  mixing  region  was  the  width 
of  the  base  region  rather  than  the  lasant  nozzle  exit  dimension. 

(2)  Comparisons  between  theory  and  zero  power  gain  and  chemiluminescence  data  showed 
that  the  model  accurately  predicted  the  experimentally  observed  trends  and  was  in  good  agreement 
with  the  absolute  magnitude  of  the  measured  zero  power  gains. 

(3)  Theoretical  comparisons  with  closed  cavity  power  data  showed  the  theory  to  overesti¬ 
mate  power  by  approximately  15%  for  most  DF  conditions  and  25%  for  HF  tests.  The  assumption 
of  plane  parallel  mirrors  (Fabry-Perot  model)  will  in  part  account  for  the  overly  optimistic  predic¬ 
tions  as  will  the  assumption  of  rotational  equilibrium. 

(4)  The  C/B  model  accurately  predicted  the  observed  variations  in  laser  performance  with 
the  dependent  flow  parameters  (mass  throughput,  temperature,  diluent  level,  and  cavity  mixture 
ratio). 

(5)  The  current  version  of  the  C/B  computer  code  has  been  shown  to  give  accurate  pre¬ 
dictions  for  low  pressure,  high  base  nozzles  for  a  wide  range  of  flow  conditions.  The  effect  of  F-atom 
wall  recombination  has  been  shown  to  be  an  important  mechanism  in  that  both  the  BCL-1 4  power/x 
and  the  LRCL-2  gain  data  are  both  accurately  predicted  assuming  complete  F-atom  recombination. 
The  results  of  the  kinetic  model  sensitivity  study  show  the  predictions  to  be  insensitive  to  factors 

of  two  changes  in  various  rates. 

(6)  The  simplicity  of  the  model  has  required  that  a  number  of  assumptions  be  made  to 
keep  the  computer  solution  time  reasonable.  Physical  processes  which  have  not  been  considered 
explicitly  but  which  have  been  or  are  currently  being  incorporated  in  various  computer  codes  at 
Bell  Aerospace  Textron  are:  (i)  rotational  nonequilibirium  effects,  (ii)  upstream-downstream 
coupling  influence  of  stable  resonators,  (iii)  transverse  pressure  gradient  effects,  (iv)  fully  viscous 
primary  and  secondary  nozzle  flow  influences  on  cavity  entrance  conditions.  Results  from  these 

studies  will  be  described  in  forthcoming  reports. 
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Figure  1.  Flow  System  and  Nozzle  Geometry 
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Figure  3.  Schematics  of  Two  High  Base  Axisymmetric  Nozzles,  BCL-13  and  BCL-14 
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Figure  5.  Sideview  of  2D  or  Axisymmetric  High  Base  Nozzle 
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Figure  7.  Variation  of  py  with  m/A 
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Figure  8.  BCL-14  -  DF*  -  Pi(J)  Lines,  Power  On 
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Figure  14.  LRCL-2,  Condition  501,  m/A  =  0.042  gm/sec-cm^ 
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o  DATA:  P2(5) 
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Condition  500,  m/A  =  0.042  gm/sec-cm^ 


Figure  17.  LRCL-1,  m  lA  =  0.147 
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Figure  18.  LRCL  Free  Jet  Performance 
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efficiencies  r,  s, . . .  respectively. 
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TABLE  3a.  REACTION  RATE^^)  DATA  FOR  THE  H,2-F2  CHEMICAL  SYSTEM  (1978); 
REACTIONS  INVOLVING  GENERALIZED  COLLISION  PARTNERS 


Reactions 

M<2) 

A 

N 

E  (kcal) 

Comments 

H+H  =  HjlV) 

*/H,  H2 

6.2  X  10” 

0.0 

V<3 

H+H  =  HjlV) 

H2,  2.0  H 

9.4  xIO"^ 

0.0 

V<3 

F+F  =  F2 

*/2.4  F2,  2.4  F 

4.71  X  10“ 

■1.25 

H2(V)  =  H2(V-1) 

*/H2,  H 

2.5x10^V 

0.0 

V<3 

H2(V)  =  H2(V-1) 

H2,H 

1.0x10  *  V 

0.0 

V<3 

H2(V)  =  H2(V-1) 

H 

2.0x10“ 

2.72 

V<3 

HF(V)  =  HF{V-1) 

HF+0.5DF 

2.95  X  10“  V 

-1.0 

0.0 

V<7 

HF|V)  =  HF(V-1) 

HF+0.5DF 

3.5x10^V 

2.26 

0.0 

V<7 

HF(V)  =  HF(V-2) 

HF+0.5DF 

7.5x  10“{V-1) 

-1.0 

0.0 

V<7 

HF(V)  =  HF(V-3) 
HF(V)  =  HF{V') 

HF+0.5  DF 

5.0x10“(V-1) 

-1.0 

0.0 

V<7 

H 

4.5  X 10“  g(v) 

0.7 

v'<,g(i)  =  0.i 
g(V>l)  =  1.0 

HF(V)  =  HF(V-1) 

Ar,  F2,  2(He,N2,CF4) 

7.7x10’’ V 

5.0 

0.0 

V<7 

HF(V)  =  HF(V-I) 

H2 

6.0x10’ V 

1.0 

0.0 

V<7 

HF{V)  =  HF(V-I) 

F 

1.93  X  10“  g(V)' 

-0.75 

3.6 

g(V)  =  V'-*’ 

HF(V)  =  HF(V-2) 

F 

9.64x10“  g(V) 

-0.75 

3.6 

g(V)=V-1)‘** 

HF(V)  =  HF(V-3) 

F 

1.45  X  10“  g(V) 

3.6 

g(V)  =  V-2)°-» 

(1)  See  Table  2a  Footnote  1 

(2)  See  Table  2a  Footnote  2 
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TABLE  3b.  EXCHANGE  PUMPING  AND  REACTION  RATE^^)  DATA 
FOR  H2-F2  SYSTEMS  (1978) 


Reactions 

A 

N 

E  (Kcal) 

Comments 

F  +  Hj  (V')  =  HF(1)  +  H 

2.72  X  10*^ 

n 

V'<3 

F  +  H2(V')  =  HF(2)-H 

8.79  X  10*® 

mim 

V'<3 

F  +  H2(V')  =  HF(3)  +  H 

4.48  X  10*® 

1.6 

V'<3 

H  +  HF(4)  =  H2  (V)  +  F 

3.70  X  10*® 

0.46 

V<1 

H  +  HF(5)  =  H2(V)+F 

3.87  X  10*®g  (V) 

0.51 

g(0)  =  1,0(2)  =  1.78 

H  +  HF(6)  =  H2(V)  +  F 

4.17  X  10*® 

0.58 

V<3 

H  +  F2  =  HF|3) 

1.2  X  10**g(V) 

mm 

2.4 

g(V)=.08,.13,  .35,  .44, 
for  V=  3, 4, 5, 6 

HF(V)  +  HF(V')  =  HF(V-1)  +  HF(V +t)  3.0x  10® 

■1.0 

0.0 

AV  =  V-V'<3 

H2  (V)  +  HF(0)  =  H2  (V-1)+  HF(1) 

8.3x  10“g(V) 

0.0 

0.0 

g(V)=  1,3.3,  10,23,46, 
90,forV  =  0, 1,2, 3, 4,  5 

HjIV)  +  H2  (V')  =  Hj  (V-1)  +  H2  (V'+l) 

3.5  X  10®  V{v'+1) 

1.5 

0.0 

V,  V'<3 

1.  See  Table  2a  Footnote  1. 

2.  See  Table  2a  Footnote  2. 
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TABLE  4.  COMPARISON  BETWEEN  THEORY  AND  DATA  FOR  THE  BCL-13 
REGENERATIVELY  COOLED  AXISYMMETRIC  NOZZLE 

(a)  Test  Conditions 


m/A,  gm/sec-cm^ 

0.249 

2721 

0.80 

0.50 

«L 

13.0 

Theory  and  Data  Comparison 

Predicted 

Measured 

Lasing  Zone  Length,  cm 

1.10 

0.86 

Relative  Power 

1.0 

0.98 
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TABLE  S.  COMPARISON  OF  MEASURED  AND  PREDICTED  X^.  FOR  LRCL  NOZZLES 


Nozzle 

m/A 

gm/$ec-cm^ 

cm 

(Xj)  C/B 

Data 

C/B 

(Xp)  Data 

LRCL-1 

0.0412 

1.77 

1.15 

0.65 

0.0988 

2.50 

0.1472 

2.45 

1.24 

0.2865 

2.98 

1.50 

0.3468 

2.95 

1.85 

0.63 

LRCL-2B 

0.0430 

0.99 

1.40 

1.41 

0.1033 

2.92 

2.00 

0.68 

0.1546 

2.94 

1.60 

0.1954 

3.41 

1.50 

0.44 

LRCL-2*’’ 

0.0405 

1.50 

0.63 

0.1002 

1.48 

0.1480 

1.99 

1.07 

0.54 

0.1920 

1.97 

1.05 

0.53 

0.2529 

2.46 

0.88 

0.36 

LRCL-2 

0.0405 

1.50 

0.80 

Q 1 

1.48 

1.70 

1.17 

0.2529 

2.46 

0.65 

(1) 

Constant  Area  Cavity  Boundary  Condition 

( V" 

(2) 

p.Pi  / ) 

Pe  Pi 

~  I  Cavity  Boundary  Condition 
0  / 
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TABLE  6.  HF  MODELING  SENSITIVITY  TO  REACTION  RATE  VARIATIONS. 
TEST  CASE  IS  FOR  LRCL-2  NOZZLE,  RUN  NO.  1  +  B6-2868 


n 

Rate  Constant  Examined 
and  How  Varied 

At  Cutoff 

o 

2Xc 

a 

1 

Baseline 

1.80 

512 

1.00 

2a 

Pumping  Reaction  Rates  x  2 

1.60 

515 

1.01 

2b 

Pumping  Reaction  Rates  x  1/2 

1.80 

497 

0.97 

3 

HF/HF  V-T  *2 

1.60 

486 

0.95 

4 

HF/DF  V-T  *2 

1.60 

494 

0.96 

5 

HF/H  V-T  *2 

1.80 

479 

0.94 

6 

HF/HF  V-V  *2 

1.60 

500 

0.98 

7 

HF/H2  V-V  *2 

1.40 

473 

0.78 

0.92 

Nomenclature: 


Barred  quantities  are  non-dimensionalized  by  their  baseline  values. 

Cutoff  is  point  where  first  four  levels  reach  maximum  power. 

(lii/A  =  0.0405  gm/sec-cm^,  Tq(.=  2126°K,  \p^  =  3.96,  Rl=  1 1.5, /Jg  =  0.0) 
EXPERIMENTAL  X^.  =  1 .50  CM,  a  =  474  kJ/kgm 
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